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Abstract 



By computing 254 unstable stationary solutions of the Kuramoto-Sivashinsky 
equation in the extensive chaos regime (Lyapunov fractal dimension D = 
8.8), we find that 30% satisfy the symmetry of the time-average pattern of 
the spatiotemporal chaos. Using a symmetry pruning of unstable stationary 
solutions, the escape-time weighting average converges to the time-average 
pattern of the chaotic attractor as 0(1/N), where N is the total number of 
unstable stationary solutions and unstable periodic orbits in the average. 
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A topic of great current interest in spatiotemporal chaos theory is to determine the 
appropriate quantities to characterize and distinguish chaotic states from one another. At- 
tempts to distinguish states using dimension densities |T|-[|, correlation lengths and 
pattern dependent order-parameters |7],||] reduce the dynamics to scalar quantities with 
varying degrees of sensitivity to parameter changes and often different parametric depen- 
dencies. Various projections of the chaotic dynamics have proven more sensitive to parameter 
changes or inhomogeneities in the dynamics. These have included histograms of the inten- 
sity of multi-mode lasers |§ , spatial averages of the curvature of a fluid surface in Faraday 
wave experiments [110], and averages of shadowgraph images in convection experiments [111 



Recent attempts to predict the structure in the projections using unstable periodic orbits 
(UPOs) in low- dimensional chaos [12|] and spatiotemporal chaos [|l3l have demonstrated that 
the escape-time weighting can approximate averages of chaos based on the moments of a 
computed set of UPOs. This suggests that UPOs may be more appropriate to character- 
ize chaos because an infinity of order-parameters (including those discussed above) can be 
approximated using a small set of UPOs. Also, the parametric dependence of the order- 
parameters are easily determined by using Newton continuation methods to evaluate the set 
of UPOs in nearby parameter regimes. 

The utility of UPOs to characterize chaos is based on the concept that Axiom-A chaotic 
systems have a dense set of UPOs associated with the strange attractor [14]. In Axiom-A 
chaotic attractors with a symbolic dynamics, trace-formulas of a small number of ordered 
UPOs can be used to estimate dynamical invariants and averages of chaos |15|,|16j. Un- 



fortunately, most spatiotemporal chaotic systems are not Axiom-A [£L7j] and do not have 
a symbolic dynamics which make trace formulas inappropriate for high- dimensional chaos. 
However, the escape-time weighting (ETW) of UPOs has been used to approximate aver- 
ages in high- dimensional chaos [^] where trace-formulas often fail primarily due to a loss 
of symbolic dynamics [Q. To date, the ETW average of the easier to compute unstable 
stationary solutions of the spatiotemporal chaos has not been evaluated. 

Our goal in this report is to demonstrate that the convergence of the ETW average of 
UPOs in the spatiotemporal chaotic Kuramoto-Sivashinsky equation can be improved using 
unstable stationary solutions (USS). We demonstrate that ETW converges as 0(1/N) with 
N computed UPOs and USS. As no method exists to enumerate UPOs and USS in nonlinear 
partial differential equations, we use utilize damped-Newton algorithms |I3] to compute the 
UPOs and USS. 

The unstable stationary solutions computed using damped-Newton methods can be iso- 
lated from the chaotic attractor. Past work computing USS in the Kuramoto-Sivashinsky 
equation with periodic boundary conditions demonstrated that the USS consisted of an un- 
expectedly large class of stationary solutions, which include laminar states, N-cell states 
(periodic patterns), long- wave modulated N-cell states, giant states (large amplitude solu- 
tions), and traveling wave states [p~8| — pO|l . Many of these stationary states were not observed 



in the chaotic dynamics of the Kuramoto-Sivashinsky equation, whereas the time- dependent 
T^,|?T| did qualitatively resemble the chaotic dynamics. 



UPOs 



Motivated by cycle-expansions [|T5|, [16 



researchers became interested in approximat- 
ing moments of sustained chaos in the Kuramoto-Sivashinsky equation using UPOs on 
the chaotic attractor [|13|,|21|]. Christiansen, Cvitanovic, and Putkaradze considered the 
Kuramoto-Sivashinsky equation with periodic boundary conditions [BTJ. Cycle-expansions 
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of the UPOs were evaluated as the small fractal dimension D < 2.1 of the system allowed an 
understanding of the symbolic dynamics of the UPOs. We consider extensive chaos in the 
Kuramoto-Sivashinsky equation (fractal-dimension D = 8.8) with rigid boundary conditions 
where the lack of a symbolic dynamics makes cycle expansions inapplicable. However, the 
escape-time weighting of UPOs |H| can be used to approximate the time-average pattern of 



the spatiotemporal chaos. 

To explore how a suitably selected set of USS improves the escape-time weighting av- 
erage of UPOs, we compute the stationary solutions of the Id Kuramoto-Sivashinsky (KS) 
equation, 

d t u = —ud x u — d\u — d^u, xe[0,L], (1) 
with 'rigid" boundary conditions u = d x u = at x = and x = L. To compute stationary 



solutions, a damped-Newton algorithm (discussed in Ref is used to find the zeros of 
the residual, 

F = ud x u + d%u + %;u, (2) 

which was discretized with Ax = 0.25. The initial guesses to the damped-Newton algorithm 
were chosen randomly from the chaotic attractor and the convergence criteria for finding 
a USS was that \\8u\\ < 10~7||it o ||, where 5u is the correction vector, u Q is the solution, 
and || • ■ ■ || denotes the infinity norm. For L = 50 and with 30,000 initial guesses, we found 
254 distinct unstable stationary solutions and no stable solutions. 

Fig. |T] shows three of the 254 USS computed by the damped-Newton algorithm. Many 
unstable stationary solutions were found to be isolated (Fig. |l](a) and Fig. |l|(b)) from the 
chaotic attractor not resembling any chaotic solution. Other USS did resemble chaotic 
solutions (Fig. [l|(c) ) and exhibited small values of the difference function, 

G(t) = \\u(t,x) - s(x)\\, (3) 

where u(t, x) denotes the spatiotemporal chaotic field, s(x) denotes the stationary solution, 
and || • • • || denotes the infinity norm. Unfortunately finding the minimum value of G(t) 
depends on recurrence times of the spatiotemporal chaos, so we prune the set of computed 
USS by requiring that stationary solutions have the same symmetry as the time-average 
pattern of the chaotic attractor. We impose the constraint that the stationary solutions 
satisfy the condition that d x u(0) < and d x u(L) > 0. This symmetry pruning results in a 
set of 80 symmetry related USS. We note that this criteria does not remove all isolated USS 
(Fig. |(a)). 



All of the average patterns of the UPOs computed for the KS equation in Ref fl3 
satisfy the symmetry pruning defined above for the USS. These average patterns were used 
to approximate the time-average pattern of the chaotic attractor using the escape-time 
weighting. The escape-time weighting (ETW), 

< m >= (4) 

approximates an average quantity < m > of the chaotic attractor in terms of the same 
average quantity rrij for each of the j computed UPOs, where r = 1/% an d \ are the 
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nonnegative Lyapunov exponents of the UPOs. The ETW of 127 UPO average patterns < 
m > qualitatively resembles the time-average pattern of the attractor < u >t (Fig. @(a)). To 
include symmetry pruned USS into the ETW average, the weighting becomes t = Ei V s » 
where Sj are the nonnegative growth factors of the stationary solution. Applying ETW to 
the spatial profiles of symmetry pruned USS and the UPO average patterns, we obtain the 
improved ETW approximation < m > to the average pattern < u >t in Fig. |2|(b). This 
improvement is quantified using the error function E(x), 

E(x) = || < m(x) > —ut(x)\\, (5) 

where || ■ • • || is the absolute difference. The average error < E(x) > for the ETW average of 
both UPOs and USS (Fig. 0(b)) is 50% smaller than that of a ETW average of just USS or 
a ETW average of just UPOs (Fig. fj(a)). Empirically, ETW appears to converge as 0(1/ N) 
with N computed UPOs and USS in a high fractal dimension chaotic system in agreement 



with the low- dimensional convergence rate of ETW reported in Ref |L2 . 

In conclusion, of 254 USS computed with a damped-Newton algorithm approximately 
30% had the same symmetry as the time-average pattern of the spatiotemporal chaos. In- 
cluding the symmetry pruned USS into the ETW average of the UPOs resulted in an im- 
provement to the average consistent with a 0(1/N) convergence of ETW with N computed 
UPOs and USS. As ETW does not require the assumptions of Axiom-A dynamics or a sym- 
bolic dynamics, it constitutes one of the few averaging techniques for high fractal dimension 
chaos. With improved numerical techniques and computational resources to compute many 
UPOs and USS, ETW may provide an important means to quickly evaluate the parametric 
dependencies of a variety of projections of spatiotemporal chaos. 

The author thanks Henry S. Greenside and Joshua Socolar for many helpful conversa- 
tions. This work was supported in part by NSF grants NSF-DMS-93-07893 and NSF-CDA- 
92123483-04, and by DOE grant DOE-DE-FG05-94ER25214. 
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FIGURES 



FIG. 1. Three representative unstable stationary solutions of the Kuramoto-Sivashinsky equa- 
tion with L = 50.0, Ax = 0.25, and rigid boundary conditions u(0) = u(L) = u x (0) = u x (L) = 0. 
The criterion for convergence was \\5u\\ < 10 _ 7||« o ||, where 5u is the correction vector, u Q is the 
solution, and || • • • || denotes the infinity norm. (A) and (B) are isolated stationary solutions, 
whereas the stationary solution (C) has finite measure on the chaotic attractor. 

FIG. 2. (A):Time-average of the Kuramoto-Sivashinsky chaotic trajectory (u(t, x))t (solid line) 
and the escape-time weighting (M) (dashed line) of 127 unstable periodic orbit time-average pat- 
terns. (B): Time-average of the Kuramoto-Sivashinsky chaotic trajectory (u(t, x))t (solid line) and 
the escape-time weighting (M) (dashed line) of 127 unstable periodic orbit time-average patterns 
and 80 symmetry pruned unstable stationary solutions. 

FIG. 3. (A): Error E{x) between the time-average of the Kuramoto-Sivashinsky chaotic trajec- 
tory and the escape-time weighting of 127 unstable periodic orbit time-average patterns of the Ku- 
ramoto-Sivashinsky equation. The average error is 0.11. (B):Error E{x) between the time-averaged 
pattern of the Kuramoto-Sivashinsky chaotic trajectory and the escape-time weighting of 80 sym- 
metry pruned unstable stationary solutions and 127 unstable periodic orbit time-average patterns 
of the Kuramoto-Sivashinsky equation. The average error is 0.05. 
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